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In this paper, we explore the possibihty of spontaneous CP violation in the scattering of quarks and 
anti-quarks from QCD Z (3) domain walls. The CP violation here arises from the nontrivial profile 
of the background gauge field (^o) between different Z (3) vacua. We calculate the spatial variation 
of Ao across the Z(3) interface from the profile of the Polyakov loop L{x) for the Z(3) interface and 
calculate the reflection of quarks and antiquarks using the Dirac equation. This spontaneous CP 
violation has interesting consequences for the relativistic heavy-ion collision experiments, such as 
baryon enhancement at high Pt- It also acts as a source of additional 3/^ suppression. We also 
discuss its implications for the early universe. 

PACS numbers: 25.75.-q, 12.38.Mh, 11.27.-|-d 



I. INTRODUCTION 

' The possibility of extended topological objects in the quark-gluon plasma (QGP) phase, e.g. Z(3) interfaces arising 
^4 from spontaneous breaking of Z(3) symmetry, has been extensively discussed in the literature [JQ- It has also been 
Q_| pointed out that there are also topological string defects in QGP forming at the junctions of Z(3) walls [!]• Formation 
D , and evolution of these objects in the initial transition to the QGP phase has been studied in the context of relativistic 
i-S^ ■ heavy-ion collision experiments (RHICE) jal . Certain consequences of Z(3) walls for baryon inhomogcncity generation 
in the universe have also been explored [6[. Investigation of these objects is important not only for probing the 
very rich vacuum structure of the QCD in the deconfining phase, but also because these provide the only example 
of topological defects in a relativistic quantum field theory which can be probed in laboratory conditions, namely, 
the relativistic heavy- ion collision experiments (RHICE). The existence of these objects has been questioned in the 
literature, especially in the presence of quarks [3, [I]- However, there are recent Lattice studies by Deka et al. @ of 
QCD with quarks which have attempted to directly probe the existence of different Z(3) vacua. These results show 
strong possibility of the existence of non-trivial, metastable, Z(3) vacua for high temperatures. The exact value of the 
temperature, above which these metastable Z(3) vacua are seen, is not important. What is important is that these 
vacua seem to exist as metastable thermodynamic phases of QCD in the deconfining regime, and hence associated 
topological objects will necessarily arise in any realistic phase transition from the confining phase to the QGP phase. 

In this paper we will investigate an interesting possibility arising from the existence of Z(3) interfaces. We will study 
refiection of quarks and antiquarks from Z(3) walls and show the existence of CP violation arising from the Z(3) walls. 
This CP violation is spontaneous, arising due to the background configuration of the gauge field corresponding to the 
Z(3) wall, and was first demonstrated by Altes et al. [lOl |. It was shown in rcf. in the context of the universe, 
that due to the non-trivial background field configuration for the standard model gauge fields, the localization of 
quarks and antiquarks on the wall is different. Its possible effects on the electroweak baryogenesis via sphalerons 
was discussed in Same possibility of spontaneous CP violation for the case of QCD was also discussed in [Tlj . 
We extend these studies by calculating the propagation of quarks and antiquarks across the Z(3) walls and show 
that they have different reflection coefficients. For this we calculate the profile of the order parameter L (x) between 
different Z(3) vacua [1| using the the effective potential for the Polyakov loop, as proposed by Pisarski [l2|- We then 
obtain obtain the profile of the background gauge field Aq from this L (x) profile. This Aq configuration provides a 
potential for the propagation of quark causing non-trivial refiection of quarks from the wall. It is important to know 
the uncertainties in the determination of the Aq profile depending on the choice of the specific form of the effective 
potential, such as those given in [l^ . To address this issue we repeat the above calculation for another choice of 
effective potential of the Polyakov loop as provided by Fukushima [l^. We find that, even though the two effective 



'Electronic address; |atreya@iopb.res.in| 
t Electronic address: ajit@iopb.res.in' 
t Electronic address: anjishnu@prl.res.in| 



2 



potentials (in refs. [1^ and [l3|) are of qualitatively different shapes, the resulting wall profile and profile are 
surprisingly similar. This gives us confidence in the use of our procedure to calculate the reflection of quark and 
antiquarks from the Z(3) interfaces. 

Different values of the reflection coefficients of quarks and antiquarks from the Z(3) walls will have very interesting 
implications for the case of RHICE and for the early universe. Here we mention that in the earlier studies by some 
of us the reflection of quarks/antiquarks from Z(3) walls (in the context of RHICE and the universe) 0], was studied 
by modeling the dependence of effective quark mass on the magnitude of the Polyakov loop, and no possibility of 
spontaneous CP violation was explored. This CP violation, resulting in different reflection coefficients of quarks and 
antiquarks from Z(3) walls, will lead to segregation of quarks and anti-quarks due to motion (collapse) of walls. As 
a result there will be selective concentration of baryon (or antibaryon) number in different regions, depending on the 
Z(3) vacua involved. This will have direct observable consequences for the relativistic heavy ion collision experiments. 
For example, it will affect the yield of baryons and mesons, enhancing baryon multiplicities and suppressing meson 
multiplicities. As we will see, these effects are expected to be important for heavy quarks, especially for charm and 
heavier flavors. A detailed analysis of these effects is planned for a future work. This CP violation can also play 
an important role in the context of early universe, especially for generation of baryon density inhomogeneities, by 
segregating baryons and antibaryons. We mention here that our analysis of reflection of quarks in this paper utilizes 
Z(3) wall profile of pure SU(3) gauge theory, without dynamical quarks. The effects of quarks may not be important 
in the context of RHICE due to small length and time scales involved, but for the case of universe these effects will 
be of crucial importance. We will discuss this further below. 

The paper is organized in the following manner. In section II we discus the basic physics of the origin of spontaneous 
CP violation due to the presence of Z(3) interfaces [13: El and discuss the effective potential for the Polyakov loop, as 
proposed by Pisarski [l^l for calculating various quantities. In section HI, we discuss how to obtain theprofile of the 
background gauge field Aq from the profile of the order parameter L (x) between different Z(3) vacua In section 
IV we address the issue of uncertainties in the determination of the Aq profile depending on the choice of the specific 
form of the effective potential by repeating the calculations of section III for the effective potential of the Polyakov 
loop provided by Fukushima jl3| . The resulting wall profile and Aq profile are found to be very close to those found 
in section III. We use the profile of Aq as calculated in section HI, for the Dirac equation (in the Minkowski space) in 
section V to calculate the reflection and transmission coefficients for quarks and antiquarks. Section VI presents our 
results and conclusions are discussed in section VII. 

II. ORIGIN OF SPONTANEOUS CP VIOLATION 

We first discuss the basic physics of the origin of the spontaneous CP violation from the existence of Z(3) walls. 
For the case of pure SU{N) gauge theory, we start with the deflnition of the Polyakov loop, [l5l - [l7j 



L(x) = —Tr 



(1) 



Pcxp^ig j Ao{x,T)dT 



where, Aq{x,t) = Aq{x,t)T"- , {a = 1,...N) are the gauge flelds and T" are the generators of SU (N) in the 
fundamental representation. P denotes the path ordering in the Euclidean time r, and g is the gauge coupling. Under 
global Z{N) symmetry transformation, the Polyakov Loop transforms as 

L{x) — > Z y. L{x), where Z = e"^. (2) 

Here, = 2Trm/N; m = 0, 1 . . . (iV - I). 

Thermal average of the Polyakov loop, {L{x)), is the order parameter for the conflnement-deconfinement phase 
transition. (From now onwards, we will use L[x) to denote {L{x)),) It is related to the free energy of a test quark in 
a pure gluonic medium {L{x) oc e~^^). L{x) ^ implies finite free energy of a test quark and hence, the deconfined 
phase (i.e the system is above the critical temperature Tc). This leads to spontaneous breaking of Z{N) symmetry. 
On the other hand, L{x) ~ implies infinite free energy of a test quark and hence, confined phase (i.e. system is 
below Tc). The Z{N) symmetry is then restored. The iV-fold degeneracy of the ground state implies the existence of 
interfaces between regions of different Z{3) vacua. For QCD, the gauge group is the color group SU{3)c- It has three 
Z(3) vacua resulting from the spontaneous breaking of Z{3) symmetry in the high temperature (deconfined) phase 
characterized by, 

L(.f) = l,e'27r/3^g»47r/3_ (3) 
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As we mentioned above, there have been questions whether these Z(3) domains have some physical meaning or not 
0, 111- The inclusion of quarks raises further issues as they do not respect the Z{N) symmetry. It has been argued 
that it is possible to interpret the effect of addition of quarks as the explicit breaking of Z{N) symmetry and lifting of 
degeneracy of the vacuum [l2l. [l8l - l20| . and we will follow this approach. Further, as we mentioned in the Introduction, 
recent lattice QCD studies with quarks @ have strengthened the physical basis for the existence of these different 
Z(3) vacua. The metastability of non-trivial Z{3) vacua will have important implications for RHICE and the early 
universe. However, for the rest of the paper we will consider the pure gauge case for calculating the Z{^) interface 
profiles. This is because our main objective here is to show the interesting possibility of spontaneous CP violation in 
the reflection of quarks and antiquarks from Z(3) walls which is independent of the explicit symmetry breaking. We 
will briefly comment on the effects of quarks in the last section, detailed study of these effects will be presented in a 
future work. 

As mentioned earlier, different Z{3) vacua have interpolating L{x) profile leading to Z(3) interfaces. This essentially 
means that there is a background gauge field Ao{x) profile which interpolates between different Z{3) vacua. The 
quarks/anti-quarks moving across the Z{3) domain walls will behave differently in the presence of a given spatially 
varying Aa field configuration. As a result, we should have different reflection and transmission coefficient for quarks 
and anti-quarks. This is the source of CP violation. The origin of this CP asymmetry is spontaneous in nature. The 
earlier studies [13, [lli of this spontaneous CP violation arising from Z(3) walls focused on the localized solution of 
Dirac equation (in Euclidean space), and it was shown that if a wave function for a fermion species localizes, then it's 
CP conjugate doesn't. The whole discussion in ref. El was within the Euclidean formalism and the exact gauge 
field profile was not determined in these investigations. 

In this paper we are interested in the calculation of reflection and transmission coefficient of quarks and anti-quarks 
and hence, in the propagating solutions. It is important to note here that the background gauge field profile comes 
from the finite temperature field theory, which is formulated in the Euclidean space. To calculate the reflection and 
the transmission coefficients (or to study propagation of quarks, in general), we need to solve Dirac equation in the 
Minkowski space. 

We start with the Dirac equation in the Euclidean space, with the spatial dependence of Aq calculated from Z(3) 
wall profile as mentioned above. Then we do the analytic continuation of the full equation to the Minkowski space 
and use the resulting equation to calculate the reflection and transmission coefficients. We should mention here that 
it may seem puzzling that we are extracting information about colored objects (i.e. ^o) starting with a colorless 
object, the Polyakov loop. However, as we will explain later in Section V, starting with a given profile of L{x), one 
does not get unique solution for Aq(x) and the ambiguity about color information manifests itself in the form of a set 
of solutions of Aq. 

We will use the effective model for the Polyakov loop as proposed by Pisarski [l2| . The Lagrangian density has the 
form 

C = ^\d,L\^T^-V{L). (4) 

= 3 for our case (i.e QCD). is mifltiplicd with the first term to give the correct dimensions to the kinetic term. 
V{L) is the potential term that has the form 

V{L) = - + {L*r) + \{\m'y,T\ (5) 

The cubic term in L{x) in the above potential, when written in terms of L{x) = |L(a;)|e*^, gives rise to cos(30) 
term that leads to three degenerate Z(3) v acua when L{x) ^ (i.e. when T > Tc). The coefficients 62, &3 and 
64, in the potential, are fixed in ref |l8l - l20l | by comparing with lattice results for the pressure and energy density 
for pure SU(3) gauge theory [HI, [H. &2 is given by 63 = (1 - 1.11/ x) (1 + 0.265/a;)^ (1 -I- 0.300/a;)^ - 0.478, where 
X = T/Tc with Tc ~ 182 MeV. The other parameters are 63 = 2.0 and 64 = 0.6061 x 47.5/16 (the factor 47.5/16 for 
&4 is to account for the additional quark degrees of freedom compared to pure SU(3) case). With the above values, 
L (x) — > y = 63/2 + i X y^bf + 4&2 (T = 00) as T — > 00. L (x) and other quantities are normalized as follows, 

L{x) — > L {x) /y, 62 — > h/y^, &3 — > h/y, h — > hy^, (6) 

so that L (x) — >■ 1 as T — > 00. The normalized quantities are then used in eqn. (O, which is then used to calculate 
the L{x) profile using energy minimization, see ref.[3l for details. Fig.l shows the plot of |i(a;)| for the interface 
between two different vacua (in the absence of quarks all the three interfaces have same profile for |L(x)|). We 
mention that the surface tension a of the Z(3) walls was estimated in refs.[Bl ^^le above effective potential and it 
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FIG. 1: Variation of |L (if) [ between different Z(3) vacua for T = 200 MeV and T = 300 MeV respectively, as a function of z. 
Note that at higher temperature, the wall thickness is smaller, as expected. 



was found that a = 0.34, 2.62, and 7 GeV/fm^ for T = 200, 300, and 400 MeV respectively. There have been Lattice 
studies of Z(3) waU tension. In ref.^l] the surface tension was found to be a{T^ = O.ITT^. With = 182 MeV the 
T = 200 resuh for a in ref.T 
for larger temperatures, T 



] is larger by almost factor 10 than the lattice result of ref.[23[. However, the values of a 
300 and 400 MeV arc in reasonable agreement with the analytical estimates [l^l (which 
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for large temperatures) 



give a 

The energy minimization program gives the full profile for L{x) which is then used for calculating Aq(x) as described 
in the next section. (As we mentioned in the Introduction, we will also consider another form of effective potential as 
provided by Fukushima [l^ in section IV.) 



III. OBTAINING Ao PROFILE 

In this section we calculate the profile form L{x) profile by inverting eqn.([T]). As in ref. [l^I we choose Aq to be 
of the form 



2ttT 

Ao = (aX-i + bXs), 

9 



(7) 



where, g is the coupling constant and T is the temperature, while A3 and As are the diagonal Gell-Mann matrices. 
Coefficients a and b depend only on spatial coordinates. The advantage of taking this gauge choice is that we are 
dealing with the eigenvalues of the matrices that are invariant under gauge transformation. 

We take Aq to be independent of t. This is for simplicity. Further, it can be justified in the high temperature limit 
due to periodic boundary conditions on Aq in the (Euclidean) time direction in the imaginary time formalism being 
used here for finite temperature field theory. 

Substituting eqn.([7]) in eqn. ([T]), we get 



3L{x) ~ exp(iQ!) + cxp(i/3) + exp(i7). 



(8) 



where, a = 27r (| + |) , /? = 27r (| — |) and 7 — 27r(— |^). On comparing the real and imaginary part of eqn. ([8]), 
we get 



cos (a) + cos (/3) + cos (7) = 3|L| cos {9) , 
sin (a) + sin (/3) + sin (7) = 3|L| sin (6*) . 



(9a) 
(9b) 



Here 9 is defined by writing L{x) = \L{x)\e^^ . In eqn. ([T]), Aq appears in the phase, so any increment in the phase 
by a factor of type 27rn will result in the same value of L {x). We first consider the above equations for L = 1 vacuum. 
Note that \L\ < 1 for finite temperatures. However, we will keep referring to the three Z(3) vacua as L = 1,Z,Z^. 
The solutions are a set of ordered pairs (a, b)j^^^. These different solution sets reflect 27rn ambiguity in Aq. Similarly, 
we find the solution sets (a, corresponding to the L = Z = exp{i2'K/i) vacuum. One now needs to find the 

appropriate values of (a, b) for the entire profile of L{x) interpolating between these two vacua. One ambiguity in 
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this is obvious. It may appear that any of the sets (a, b)j^^^ could be matched to any of the sets (a, h)^^^ as all sets 
for a given vacua are equivalent. However, this could lead to different Aq profiles in between, which in turn would 
lead to different reflection and transmission coefficients. This problem is resolved when we realize that the variation 
of ^0 should be smooth across the domain wall. Thus, we can simply start with any one pair (a, b)j^^^, and set it as 
the initial condition for the generation of the profile of Aq as one traverses the wall starting from L = 1 vacuum to 
L = exp{i2Tr/3) vacuum. We only require that a and b vary smoothly as the profile of L{x) changes smoothly across 
the wall. It will then automatically lead to the appropriate values of (a, h) j^^^ as L = Z vacuum is approached. 

For the results shown here we had taken the initial values of (a, b) — (—1.5, —1.0) for the L = 1 vacuum (in a region 
far left to the interface). As one approaches the interface, say, along the z axis, new value of L{x) is selected from 
the profile of L{x) (calculated from the energy minimization program). We then take small range of values near the 
original (a, b) = (—1.5, —1.0) and L{z) was then calculated for all these values. Those values of a and b were selected 
for which the error between the calculated L and L obtained by energy minimization was minimum. The process was 
then repeated for each value of z to obtain a, & values. Comparison between the calculated \L\ profile and the one 
obtained by energy minimization is given in figure ([2^). It clearly shows that this technique works well. Figure 
shows profile of parameters a and b across the domain wall. 
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FIG. 2: On left: Plot of calculated \L\ and the one obtained from minimizing the energy. The inset figure shows the deviation 
between the two profiles. On right: Variation of a and b between the regions L{x) = 1 and L{x) — e^^^^'"^. Initial point is 
(-1.5,-1.0) 



The calculated a, 6 were then used to calculate Aq using cqn ([7]). The Aq profile thus obtained is reasonably well 
fitted to the function Ao{x) = ptanhlqx + r) + s using gnuplot. The calculated Aq profile and fitted Aq profile are 
plotted in figure (jS]). 



IV. CALCULATION OF Aq PROFILE FOR A DIFFERENT EFFECTIVE POTENTIAL 

We now address the issue of the uncertainties in the determination of the Aq profile depending on the choice of the 
specific form of the effective potential. Other parametrization of the effective potential for the Polyakov loop have 
been given in the literature, e.g. in refs.fH, [iJI, and we will repeat the calculations of the previous section for the 
effective potential of the Polyakov loop as provided by Fukushima [l3| . For spatially varying L configurations, we will 
continue to use the derivative terms as in Eq.(4) with general dimensional considerations (with suitable normalization 
of L). The effective potential for ref.[l^ has the following form 

V[L]/T^ = ~2{d - l)e-'"''/^|TrLp - ln[~\TrL\^ + 8Re{TrLf - 18|rrip + 27] (10) 

a = (425 MeV)2 is the string tension and 2{d - \)e-'^°'l^'i = 0.5153 with Td = 270 MeV is taken as the transition 
temperature by choosing the lattice spacing a — (272 MeV)""'^. Note that for consistency with the notations of 
Rcf. [13], we will use and interchangeably, both meaning the deconfinement transition temperature. L is the 
Polyakov loop but without the normalizing factor of (= 3). (Thus, using with Eq.(4) we re- write the above effective 
potential in terms of the normalized Polyakov loop. Henceforth by L even for the above equation we will mean this 
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FIG. 3: Plot of calculated Aq and the fitted profile {Ao{x) = pta,nh{qx + r) + s). The parameters have values p = —378.27, 
q = 7.9500f, r — —49.7141, s = —1692.48. Only (f, f) component of Ao is plotted. The other components also have similar fit. 



normalized Polyakov loop). It has been argued by Schaefer et al. [25| that the transition temperature has to be tuned 
depending on the number of quark flavors Nf (and also the value of the baryon-chemical potential). In ref. [25|, the 
value of Trf = 270 MeV corresponds to the pure SU(3) case with Nj = 0. In section II we have used the effective 
potential where the coefficient 64 is suitably iiormalized for the case of 3 flavors, Nf = 3. For the case of Nf = 3, the 
value of transition temperature from ref. [23] is Td ~ 178 MeV. Thus, we will use this value of Td for the effective 
potential in Eq.(lO). 

The effective potential in Eq.(lO) is of qualitatively different nature than the one given in Eq.(5). For small values 
of L the two forms will be similar as one can see by the expansion of the Logarithmic term in the above equation. 
However, for \L\ approaching 1 the two potentials are dramatically different. V[l] in Eq.(lO) diverges at this limiting 
value thereby constraining \L\ within value 1. There is no such constraint in Eq.(5). Even the shape of V[L\ is very 
different away from the origin, especially near the three Z{3) vacua. It is thus reasonable to expect that the resulting 
profile of Z{3) wall and resulting Aq profile (using calculations of previous sections) for Eq.(lO) may be quite different 
from the ones obtained in section III for Eq.(5). 

With diverging V[L] at \L\ = 1 in Eq.(lO), and due to its non-trivial shape near the Z{3) vacua, the application 
of the technique of ref.[3l determination of L profile between two Z(3) vacua is much more complicated here. 

Especially non-trivial is the choice of initial ansatz for the wall profile which is used for the energy minimization 
program. In ref.Q, the initial profile was taken to linearly interpolate between the two Z(3) vacua as a function of 
spatial coordinate z. This choice simply does not work for Eq.(lO) due to the fact that V[L] diverges at \L\ = 1 and 
linear interpolation takes it outside this bound. For this we chose the initial trial profile to consist of two parts, one 
linearly decreasing (with z) to L = along 9 = from the vacuum value and join this with the second part linearly 
increasing (with z) along = 27r/3 to the second vacuum value. This keeps the initial profile within the allowed region 
oiV[L] in Eq.(lO). 

Second complication arises with the algorithm of energy minimization itself. In ref. Q correct L profile was obtained 
from the initial trial profile by fluctuating the value of L at each lattice point and determining the acceptable fluctuation 
which lowers the energy (with suitable overshoot criterion etc. as described in detail in ref.Q)- However, with Eq.(lO), 
fluctuations of L can take it out of the allowed region of V[L]. For this, we skip those fluctuations which take L 
outside the allowed region. With these modification in the procedure, we were able to determine the profile of the 
Z(3) wall and associated Aq profile. In section HI we had calculated the profiles for temperature T = 400 MeV (with 
Tc = 182 MeV for the effective potential in Eq.(5)). For the sake of comparison with that case, for V[L] in Eq.(lO) 
with Tc = 178 MeV [H], we calculate the profiles for T = 391 MeV which is close enough to the value T = 400 MeV, 
and has the same value for T/Tc- 

Fig. 4a shows the wall profile of \L\ for V[L] in Eq.(lO) (again, with normalized L). The profile is almost the same 
as the one shown in Fig. 2a. We mention here that for Fig. 4a we have used the same value of the coefficient of the first 
jTrLp term in Eq.(lO) as with Td = 270 MeV (by suitably changing the values of string tension etc.). This is so that 
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FIG. 4: (a) Plot of the profile of \L\ corresponding to the effective potential in Eq.(lO). (b) Comparison of the profiles of \L\ 
for different choices of Td in Eq.(lO). 




the shape of the barrier near the confining vacuum remains unaffected (which determines the first order nature of the 
transition). In any case, the overall features of the profile of the wall, such as its width and height, should depend 
more on the temperature scale rather than on the shape of the barrier for the confining vacuum. To check this, we 
also calculate the wall profile of \L\ for Eq.(lO), but now with the value of T<j = 270 MeV and T = 400 MeV. The 
comparison of the two profiles is shown in Fig. 4b. We see that the two profiles are very close to each other confirming 
above arguments. 

We recalculate the plots of a and b for the case with T = 391 MeV (with Td = 178 MeV). The resulting plots are 
shown in Fig. 5a which are seen to be very similar to those on Fig. 2b. Finally, the profile of Fig. 5b is also very 

close to the one in Fig. 3. Note that though overall all the plots in Figs. (4), (5) are very close to the corresponding 
plots in Figs. (2), (3), there is one clear difference. The profiles in Figs. (4), (5) have somewhat sharper variations from 
their asymptotic values compared to the case in Figs. (2), (3). This originates from the qualitatively different shapes of 
the two potentials in Eq.(5) and Eci.(lO) near the region of Z(3) vacua, and in that sense characterizes the difference 
in the two potentials. 

These results are quite remarkable. Even though the two effective potentials Eq.(5) and Eq.(lO) (from refs. [l^ 
and [3l) are of qualitatively different shapes, the resulting wall profile and Aq profile are almost the same. As we 
mentioned above, for small values of L the two effective potentials will have similar forms, which are fitted with the 
Lattice data. Our results thus point out that the profile of L (and consequently, the profile of A'^) are primarily 
determined by the small L region of the effective potentials. This is likely to happen if the variations near the Z(3) 
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vacua arc primarily in the magnitude of L and not in its phase. The robustness of our results against different choices 
of the effective potentials gives us confidence in the use of our procedure to calculate the reflection of quark and 
antiquarks from the Z(3) interfaces. Since the Aq profiles of Fig. 3 and Fig. 5 are almost the same, the resulting values 
of reflection coefficients for quarks/ antiquarks will also be very similar. In the rest of analysis in the paper, we will 
use the effective potential as given in Eq.(5). 



V. CALCULATING REFLECTION AND TRANSMISSION COEFFICIENTS 

To calculate the reflection and transmission coefficient, we need the solutions of Dirac equation in the Minkowski 
space. We start with the Dirac eqn. in the two dimensional Euclidean space 

[i7°5o<5^' - .97e^'o' (^) + {illdz + m)5=^] i'k = 0, (11) 

where 7° = 17° and 7^ = 7'^ are the Euclidean Dirac matrices, do denotes d/dr with t ~ it being the Euclidean time, 
j, k denote color indices. We now analytically continue the eqn (|lip to the Minkowski space to get 

[i7°ao(5J''= + 57°^'o''(-z) + («7^% + m)S''']i[^k = 0. (12) 

where now do denotes d/dt in the Minkowski space. Note that the Aq in eqn ([T2|) . which is in the Minkowski space, 
is fundamentally different from the ^0 in eqn (jlip which is in the Euclidean space. However, it is the same domain 
wall profile (i.e same Aq dependence on z) that appears in both the cases, which is what is needed for the calculation 
of reflection and transmission coefficients. For a wave function with time dependence ^p{x)e~'^^* , the eqn (fT2l) reduces 
to 

[7V535^''' + j°mP'']Mx) ^{E- Vo{z))Mx)- (13) 

where V{z) ~ —gA!i^{z) is the potential as seen by the incoming fermion. We do not have any analytic way to calculate 
the reflection and transmission coefficients for a general smooth potential, so we follow a numerical approach. Kalotas 
and Lee [26} have discussed a numerical technique to solve Schrodinger equation , approximating a general smoothly 
varying (in space) potential in terms of a sequence of step functions. We follow their approach and apply their 
technique for solving the Dirac equation (eqn (llSp ). We approximate the actual potential by n step potentials in 
series, each of equal width w as shown in figure ©. Let tpj be the wave-function for the j*^ bin and the height of 





V{z)- 
















*■ 

1 J 






1 1 \y^~^ 
1 ' 1 




1 ♦ 






(Zo + nw,0) z 



FIG. 6: Potential (V(2)) approximated by a sequence of n step potentials, each of width w. 



potential be Vy (Wc consider spin up wave function and restrict to no-spin-flip situation.) The height of the j"' step 
potential is taken to be the mean value of V{L + jw) and V(L + (j + l)w), i.e 

^ \V{L^2W)-rV{L^{2^\)w)\ ^^^^ 
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We now apply boundary conditions at j**^ step i.e at z = L + jw. This gives us a set of two equations, which when 
iteratively solved give 



M{L + jw,kq) 



Ar\L, kin) X M(L, ki)x ... M-\L + nw, k^) x M{L + nw, fcout) ' 



gikq(L+jw) 
.ifc,(i,+j™)j,^ 
E„+m 



Ea+m 



(15a) 
(15b) 



with kq = \J E"^ — iTT-^, and Eq = E — Vq. (Here no left moving wave is allowed in the region far right of the interface.) 
The reflection and transmission coefficients are then given by 



R = 
T = 
where r = 



/ret 




Bin 


Jin 




^in 



'J'm 
^out 



E + m 



E-Vm 



(16a) 
(16b) 
(16c) 



Here, kin = ^/ E"^ - m? and kout = \/{E - VoY - m^. 



VI. RESULTS 



We first calculated the reflection and transmission coefficients by assuming the Aq profile to be a step function 
rather than a smooth one, with the height of the step function being the same as that of the interface in Fig. ([6|). 
In this approximation one can calculate the reflection and transmission coefficients analytically. For anti-quarks the 
refiection and transmission coefficients are obtained by changing g — >■ —g, as anti-quarks are in 3 representation of 
SU (3). We have chosen the energies of the particles such that i? > y -I- to, so as to avoid the Klein paradox regime. 
Note that ii E < V (but V — E < m so that one is away from Klein paradox situation), then the reflection coefficient 
for quarks is 1 (repulsive potential) but for antiquarks reflection coefficient will be very small with —1/ providing 
the attractive potential. This will provide the most dramatic difference between the reflection of quarks and that of 
antiquarks from Z(3) walls. However, for the relevant energies of quarks/antiquarks at RHICE, we discuss in detail 
the case with E > V + m. 

The results for different quarks and anti-quarks (with E = 3.0 GeV for each case) are given in table ID It is clear 
that quarks have different reflection coefficients than their CP conjugates. Also, the effect is significantly higher for 
the heavier quarks (for example charm quark). 





U 


d 


s 


c 


E{GoV) 


3.0 


3.0 


3.0 


3.0 


m(MeV) 


2.5 


5.0 


100 


1270 


Rq 


1.73 X 10^'^ 


6.76 X 10~^ 


2.8 X lO"'* 


0.14 


Rq 


1.92 X 10~* 


7.55 X 10-** 


3.2 X 10"^ 


6.5 X 10"^ 



TABLE I: Table for the reflection coeflicients for various quarks in the step function approximation. Reflection is higher for 
heavier quarks. 



We now calculate the reflection coefficient for charm quark using the exact potential. The product of the matrices in 
eqn were calculated by a FORTRAN code and also by using Mathematica. Eqn (fT6|) were then used to calculate 
the reflection coefficient. At i? = 3 GeV, we get R = 0.0011 for c quark while for c the result is i? = 5.24 x lO^^*'. As 
an additional check on the results (for the smooth profile), we consider shrinking of the profile of Aq in z direction, 
and compared the reflection coefficient (for the c quark with 3 GeV energy) with the step potential result. The results 
are summarized in Table |TT1 We see that the numerical results approach the analytical results of the step function as 
Aq profile is shrunk along z to better approximate a step function. This gives us the confidence that our numerical 
technique of solving the Dirac equation is reliable. 
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Shrinking Factor 


Reflection Coeff 


No shrinking 


0.0011 


0.5 


0.017 


0.05 


0.119 


0.005 


0.123 


Step Potential 


0.140 



TABLE II: Table for the reflection coeflicients for c quark, with 3 GeV energy, when the profile is shrunk. Results approach 
the step potential as the profile gets narrower. 



It is clear that if one considers the situation of quarks/antiquarks coming from right in Fig. 1^ (i.e. approaching the 
domain wall from the side with L — Z) then antiquarks will have larger reflection coefficients while quarks will have 
smaller reflection coefficients. Also we should mention that Eqn is solved by using one component of Aq profile 
{A}f^ in this case), which gives us the reflection coefficient for one particular color (say red). Reflection coefficient for 
other colors will remain the same when SU(3)c gauge transformation is applied on the quark as well as on the vector 
potential. However, there is still an ambiguity of starting with different initial sets (a, 6) (say in the L = 1 vacuum). 
Different sets lead to different profiles for (a, b) across the domain wall, thus profile depends on the initial condition 
(which, in turn, will lead to different reflection coefficients for a quark of a given color). 

As we mentioned earlier, this ambiguity is reasonable in view of the fact that we are extracting information about 
a colored object (Aq) starting from a colorless variable L(x). Thus there is no reason to expect unique solution for 
Aq starting from a given L(x) profile, even in the diagonal gauge where is determined in terms of real (a, b). 

For several sets of values of (a, b) we have checked that different choices of (a, b) are related to each other by color 
transformation. We can explain it in the following way: Say we start with (oi, bi) for L = 1 vacuum and calculate the 
profile {a{x),b{x)) leading to profile of Aq. Now AI^,Aq'^,Aq^ all have different profiles and correspond respectively 
to scattering of red, blue, and green quarks respectively, from the given domain wall profile. Now if we start with 
a different set (02,62) and calculate the profile of Aq then we find (for example) that new ^q"'^ is the same as old 
Ag^ (where one started with (ai,6i)) and new A^^ is the same as old ^q^. This means that (02,62) set gives same 
reflection for blue quark as (oi,6i) gives for the red quark. Thus we say that our different choices of (0,6) amount 
to considering quarks of different colors for a given domain wall profile. Or, equivalently, for the scattering of a fixed 
color (say red) quark, different sets (o, b) lead to domain wall profiles carrying different color information. (We should 
mention that this holds for many sets (o, b) we have checked. However, we do not have a general proof that this 
should be true for all sets, though it looks very likely in view of the above arguments). 

For example, if we start with (o', 6') = (o, —b), i.e with (—1.5, 1), then eqn ([7]) tells us that Aq^ = and N^"^ = 
A^ . See, Fig ([7]) for the corresponding profile of (0,6). In color space Aq is diagonal with elements {A^ ,A^,A^^)^ 
and it acts on the color triplet (r, 6, g)^ . So, acting on (1, 0, 0)"^ is same as A!q^ acting on (0, 1, 0)"^ which is same 
as making different choices in color space. 

So, the ambiguity related to various (o, 6) profiles or, equivalently, corresponding ^0 profiles, seems to be the artifact 
of the ambiguity of making a color choice for the domain wall profile in terms of Aq^ starting from the domain wall 
profile in terms of L(x). 

This raises an important question whether we should be dealing with colored domain wall profile (given in terms 
of Aq profile) at all, or we should restrict to colorless objects like L(x) (which is what was done in our earlier works, 
see rcf.4,6). After all, the effective potential which we use is given in terms of L(x). Here we think that there is no 
reason to restrict to colorless objects. We are dealing with the QGP phase and there is no requirement of physical 
observables to be color singlets. If we were dealing with the confining phase then we had obligation of dealing with 
colorless objects as physical observables. For QGP phase, it should make perfect sense to think of domain wall profile 
having color properties as it is arising from Aq profile. Of course it is possible that actual domain wall profile is 
color insensitive, and quarks of all colors have same refiection coeff. from a given wall. But it is also possible that 
wall is colored and a given wall has different reflection for quarks of different colors. The only requirement of gauge 
invariance is that when color gauge transformations are done on ^0(2;) profile as well as on quarks, then numbers 
should not change, which is obviously true with the Dirac equation we are using. 

VII. DISCUSSION 

This CP violation will have interesting observable consequences for the Relativistic Heavy Ion Collision experiments 
at RHIC and at LHC. If QGP is formed in these experiments (and there are strong indications of that), then various 
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FIG. 7: On left: Variation of a for different initial values of a, b. As a is unchanged, it's profile is unaffected. On right: Variation 
of b for different initial values for b. As b changes sign in the initial values, it's profile also changes. 



Z{3) domains will inevitably be formed, leading to the formation of Z(3) walls. (Wc mention that the QGP strings 
which also necessarily form during transition to QGP phase should also lead to spontaneous CP violation. Its effects 
on quark/antiquarks scattering, or possible localization on the QGP strings needs to be explored). As these domain 
walls move/collapse, quarks/anti-quarks will get reflected/transmitted differently from these domain walls leading to 
the segregation of quarks and anti-quarks. The concentration of quarks (or antiquarks, depending on the collapsing 
vacuum) will grow in different regions of the QGP. As the effects would be stronger for heavier quarks (Table |l|, this 
should lead to enhancement of strange and charmed baryons along with the suppression in the yield of corresponding 
mesons (such as J/ip)- 

Detailed exploration of the formation and evolution of Z(3) walls and QGP strings in the context of RHICE has 
been carried out in ref. @. These simulations show that in the typical region of QGP formed in RHICE, one expects 
several Z(3) domain walls to form, their numbers ranging from 1 to 4,5. The walls may extend throughout the QGP 
region with size of order 10 fm. There are closed domain walls formed with initial size of about 5-8 fm. The velocities 
of these walls was also estimated in ref. [5| and were found to range from 0.5 to 0.8. For detailed discussion of the 
properties of Z(3) wall and QGP string networks expected in RHICE, see ref.Q. These results about the sizes and 
numbers of Z(3) walls and QGP strings are very important. This is because one should realize that in a very large 
sized QGP region, as in the early Universe, for every domain wall connecting 9 = and 6 = 2tt/3 vacua, there will 
be one connecting 6 = and = 47r/3 vacua. These walls are conjugate of each other and the reflection of a quark 
from the first wall is identical to the reflection of an antiquark from the second wall. These two walls are strictly 
degenerate, even in the presence of explicit symmetry breaking effects from dynamical quarks. Thus, on the average 
there will not be any bias for quarks and antiquarks as they scatter from a network of Z(3) walls. 

This is, however, not true for a small QGP region as produced in RHICE. As the number of Z(3) walls produced 
in such a small region is of order one there may be a net effect for the concentration of baryon number, or for 
anti-baryon, in each event. This can be revealed by event-by-event analysis. Even statistically, for a large number of 
events, one can calculate the variance of baryon number density, and spontaneous CP violation from Z(3) walls may 
be detected. For a given event also, segregation of baryons and antibaryons will occur over large distances of order 
several fm as indicated by the typical wall size and separation [5|. 

This CP violation can also be very important in the context of early universe where it can have interesting implica- 
tions for generation of baryon inhomogeneities. As collapsing domain walls preferentially sweep quarks (or antiquarks), 
segregation of quarks and antiquarks will occur. One can then discuss the formation of baryonic (or antibaryonic) 
lumps. These baryon inhomogeneities can be of large magnitude, with large separations in the context of certain low 
energy inflationary models (but now with CP violation incorporated). We will present a detailed study of this in 
a future work. 

Another important consequence will be on the Pt spectra of hadrons. The quarks/anti-quarks with high momenta 
will undergo non-trivial scattering from these Z(3) walls. As Z{3) walls collapse, some get transmitted while others 
are reflected back. For Z(S) walls forming closed, collapsing, structures, the quarks suffer multiple reflections inside 
the wall, resulting in an increment in their transverse momenta. This process continues until the walls either melt 
away or collapse completely. So the final transverse momentum of some quarks may be reasonably enhanced before 
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they escape. One can then use a specific model (such as Recombination/Coalescence model) to study the Pt spectra 
of final state hadrons, which should show an increase in the yield of hadrons at high Pf. This has been discussed 
in ref.fi^l, however, no account of CP violation was considered in that work. In the presence of CP violation, the 
modified Pt spectra will be different for quarks and for antiquarks. We plan to carry out these analyses in a future 
work. 

The most important limitation of our analysis is the absence of quark effects. Dynamical quarks will lead to lifting 
of degeneracy between different Z(3) vacua, making L = 1 vacuum as the true vacuum as discussed in refs. [T2|. [T8l - [20| . 
The one-loop corrections from dynamical quarks have also been discussed in refs.[2^-j3l|. As we mentioned, recent 
lattice studies Q have provided evidence for the existence of such metastable Z(3) vacua. Our analysis above of 
calculation of Aq profile and calculation of reflection coefficients for quarks and antiquarks can be straightforwardly 
applied for this non-degenerate case and work is underway on this. Apart from affecting the numbers (for reflection 
coefficients), its most important effect will be on the evolution of Z(3) wall and QGP string network, (see ref. [s^] for 
a detailed simulation study of these aspects). However, for the case of RHICE, due to small length (and time) scales 
involved, the dynamics of Z(3) walls is likely to remain dominated by the surface tension effects with the difference 
in pressure between different vacua not playing dominant role for such length scales). Thus the above mentioned 
features of effects on hadron spectra due to CP violation may remain qualitatively true for RHICE. 

However, for the universe the entire issue of formation and evolution of Z(3) walls crucially depends on the im- 
portance of quark effects. Some discussion of this has been provided in @ and we plan to investigate these issues in 
future in detail. Most important issue will be to see whether the spontaneous violation of CP discussed here can lead 
to a net separation of baryons and antibaryons in the universe which will have observational consequences (e.g. from 
the strongly constrained nucleosynthesis, which can be used to constrain various parameters of the model.) 
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